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Chapter  1.  Introduction 

Ferromagnetic  shape  memory  alloys  (FSMAs)  have  attracted  a  strong  interest  among 
actuator  designers  as  possible  fast-responsive  compact  airborne  actuator  material.  Earlier 
works  on  FSMAs  were  focused  on  phase  transformation  of  Fe-Pd  system  [1]  and  Ni-Mn- 
Ga  system  [2]  under  temperature  change.  There  are  three  mechanisms  of  actuation 
associated  with  FSMAs,  under  magnetic  field  which  can  be  used  as  driving  force  for  fast 
responsive  actuator  materials  [3-7]:  (i)  magnetic  field-induced  phase  transformation,  (ii) 
martensite  variant  rearrangement  and  (iii)  hybrid  mechanism. 


Figure  1.  Three  dimensional  phase  transformation  diagram  of  a  ferromagnetic  shape  memory 
alloy(FSMA)  under  applied  three  loading  of  stress,  temperature  and  magnetic  field. 

The  first  mechanism  is  based  on  the  phase  change  from  austenite  to  martensite  under 
increasing  magnetic  field,  or  reverse  phase  transformation  under  decreasing  magnetic 
field.  If  we  can  construct  three  dimensional  phase  transformation  diagram  of  stress  (c)- 
temperature  (T)-magnetic  field  (H)  axes,  see  Fig.  1,  the  phase  transformation  of  austenite 
to  martensite  phases  under  modest  magnetic  field,  requires  that  the  T-H  phase  boundary 
surface  be  inclined  toward  the  T-axis  with  the  T-H  boundary  line  that  has  a  smaller  angle 
with  T-axis,  otherwise,  the  increasing  H  loading  would  not  intersect  the  T-H  phase 
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boundary  surface.  Kato  et  al  [8]  made  preliminary  estimate  of  the  necessary  magnetic 
energy  to  induce  a  phase  transformation  based  on  thermodynamic  model,  to  conclude  that 
large  H  field  is  required  for  the  phase  change  to  take  for  both  NiMnGa  and  FePd. 
Therefore,  this  mechanism  is  not  suited  for  use  in  designing  compact  actuators  which 
may  need  a  small  and  portable  electromagnet  system  as  a  driving  unit. 

The  second  mechanism  is  to  induce  the  strain  in  a  FSMA  with  100  %  martensite 
phase  subjected  to  constant  H-field  which  acts  on  the  magnetic  moments  in  magnetic 
domains  that  exist  in  the  martensite  phase  so  as  to  rotate  them  along  the  easy  axis,  i.e.,  c- 
axis  in  the  case  of  NiMnGa  and  FePd.  The  strain  induced  by  this  mechanism  is  a 
function  of  c/a  ratio  of  FSMA,  i.e.,  the  order  of  shear  strain,  given  by  a/c  -  c/a.  Thus, 
smaller  the  c/a  ratio,  the  larger  shear  strain  can  be  induced  by  this  mechanism.  The  c/a 
ratio  of  NiMnGa,  is  reported  to  be  0.94,  which  could  provide  6  %  or  more  strain. 
Recently,  Yamamoto  et  al  the  reversible  straining  of  a  single  crystal  FePd  under  a 
constant  magnetic  field  where  the  boundary  between  two  kinds  of  variants  are  observed 
to  move,  contributing  to  reversible  strain  of  up  to  0.49%  [9].  Even  though  the  strain 
induced  by  the  second  mechanism  is  large,  the  corresponding  stress  remains  to  be  modest 
as  several  MPa  under  modest  applied  magnetic  flux  density  (IT). 

The  third  mechanism  which  we  call  as  “hybrid  mechanism”,  is  based  on  a  set  of  chain 
reactions,  first  applied  magnetic  flux(or  field)  gradient,  magnetic  force,  stress  induced 
martensite  phase  transformation  (see  Fig  1),  resulting  in  the  phase  change  from  stiff 
austenite  to  soft  martensite  phase,  leading  to  large  displacement.  The  advantages  of  this 
are  large  stress  (hundred  MPa  in  the  case  of  FePd),  modest  -  intermediate  strain,  fast 
actuation  time.  This  phase  change  can  be  applied  by  approaching  a  compact  and  portable 
magnet  close  to  the  FSMA  specimen  which  provides  a  large  magnetic  field  gradient,  thus, 
suited  for  use  in  designing  actuators  with  large  force  capability.  Currently,  the  cost  of 
processing  of  known  FSMAs  is  very  high  and  their  mechanical  behavior  is  often 
unsatisfactory  (poor  ductility).  Also,  the  processing  techniques  are  quite  cumbersome 
(making  single  crystals  FePd,  NiMnGa  etc).  To  overcome  these  difficulties,  we  proposed 
use  of  FSMA  composites  with  the  hybrid  actuation  meChanism[10]  where  two  cases  of 
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FSAM  composite  design  were  discussed,  laminated  composite  for  bending  type  actuators 
and  spring  with  complicated  cross  section  made  of  ferromagnetic  material  and  super 
elastic  SMA. 

The  advantages  of  the  proposed  FSMA  composites  are,  fast  response,  large  force  and 
stroke,  yet  its  actuators  can  be  made  into  compact  and  light  weight  with  modest  material 
and  process  costs.  As  a  demonstration,  we  designed  a  membrane  FSMA  composite 
actuator  for  use  as  a  new  synthetic  jet  actuator  which  provided  fast  jet  stream  with  speed 
of  190m/sec  [11], 

This  final  report  states  the  modeling  of  FSMA  composites  in  the  first  two  chapters 
(Chaps.  2  and  3),  followed  by  the  processing  of  laminated  and  particulate  composites  in 
chapter  4  and  5,  respectively.  Finally,  concluding  remarks  which  include  future  work  to 
be  recommended  in  Chapter  6.  Extensive  list  of  references  are  given  at  the  end,  followed 
by  the  details  of  the  analytical  modeling  in  Appendix  A,  and  list  of  publications  under 
this  project  are  given  in  Appendix  B. 
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Chapter  2.  Design  of  FSM A  Composites  by  Simple  Modeling 

In  this  modeling  study,  two  cases  of  loading,  bending  and  the  twist  modes  of  the 
composites  are  considered  with  emphasis  on  how  the  geometry  and  the  mechanical 
properties  of  the  components  influence  the  superelastic(SE)  shape  memory  alloy(SMA) 
behavior  of  the  composite.  First,  the  bending  deformation  of  the  composite  plate  with 
application  to  torque  actuators  is  theoretically  analyzed.  That  is,  the  relation  between  the 
curvature  and  the  bending  moment  for  the  composite  plate.  Next,  the  spring  of  the 
composite  wire  with  the  rectangular  section  form  is  designed  in  consideration  of 
application  to  spring  actuators,  and  the  deformation  characteristic  of  the  spring  is 
examined.  For  both  models  of  bending  and  torsion  of  FSMA  composites,  the  optimized 
microstructures  of  the  composites  are  identified. 

2.1.  Modeling  on  Superelastic  behavior  of  bending  composite  plates 

For  bending  type  actuation,  the  laminated  composite  plate  composed  of  a 
ferromagnetic  material  layer  and  superelastic  SMA  layer  as  shown  in  Figure  2(a),  is 
examined.  The  composite  plate  is  subject  to  bending  moment  M  induced  by  the  magnetic 
force  generated  by  the  ferromagnetic  material.  After  the  maximum  bending  stresses  on 
the  plate  surface  of  SMA  layer  reach  the  transformation  stress  (onset  of  superelastic 
plateau  in  the  upper  loop  of  the  stress-strain  curve,  Figure  3(b)),  the  phase  transformation 
proceeds  from  the  plate  surface  as  shown  in  Figure  2(b).  The  stress  in  the  transformed 
region  remains  constant  due  to  the  superelastic  behavior  of  SMA.  It  is  assumed 
throughout  in  this  paper  to  facilitate  the  analysis  that  the  superelastic  loop  of  SMA  is 
"flat"  i.e.  no  working-hardening  type  slope  allowed,  and  the  Young's  modulus  of  the 
austenite  is  the  same  as  that  of  the  martensite.  These  assumptions  would  allow  us  to 
obtain  closed  form  solutions  in  the  present  model,  although  the  predictions  are  still  to  the 
first  order  approximation.  The  aim  of  using  this  simple  model  is  to  identify  the  best 
thickness  ratio  of  a  ferromagnetic  layer  and  SMA  layer  in  the  composite  plate. 

Then,  the  relation  between  the  bending  moment  and  the  curvature  is  theoretically 
calculated  by  using  stress-strain  curves  of  the  constituent  materials.  Figure  3(a)  shows  the 
analytical  model.  Radius  of  curvature  of  the  composite  plate  subject  to  bending  moment 
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M  is  p,  the  thickness  of  the  composite  plate  is  h,  the  thickness  of  the  ferromagnetic  layer 
is  hf,  and  the  plate  width  is  b.  Figure  3(b)  shows  the  stress-strain  curves  of  the 
ferromagnetic  material  and  the  superelastic  SMA,  where  the  Young's  modulus  of  the 
ferromagnetic  material  is  Ef,  that  of  the  SMA  are  Esma>  the  yield  stress  of  the 
ferromagnetic  material  is  Of,  and  only  elastic  portion  of  the  ferromagnetic  material  is 
shown.  The  onset  stress  for  phase  transformation  of  superelastic  SMA  is  Go,  the  onset 
stress  for  reverse  transformation  is  Oi  in  the  superelastic  loop  portion  of  SMA.  As  a 
result,  the  relation  between  the  bending  moment  and  the  curvature  of  the  composite  plate 
also  is  expected  to  exhibit  the  superelastic  loop  if  properly  designed.  This  superelastic 
loop  of  the  ferromagnetic  shape  memory  alloy  composites  is  indeed  desired. 


(a)  Plate  bending  model  (b)  Stress-strain  relations  for  Ferromagnetic 

material  and  superelastic  SMA 

Figure  3.  Material  properties  and  model  for  the  theoretical  examination. 
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The  curvature  which  reaches  yield  stress  Of  in  a  ferromagnetic  layer  and  the  curvature 
which  reaches  transformation  stress  Oo  in  superelastic  SMA  layer  are  strongly  influenced 
by  the  mechanical  properties  and  the  thickness  of  both  materials.  Stress  distribution  is 
classified  into  the  following  three  cases  because  of  the  relation  between  the 
transformation  stress  in  the  SMA  layer  and  the  yield  stress  of  a  ferromagnetic  layer. 

Case  1 :  The  stress  in  a  ferromagnetic  layer  reaches  the  yield  stress  Of,  before  reaching  the 
transformation  stress  Oo  in  the  superelastic  SMA  layer. 

The  stress  distribution  of  this  case  upon  loading  and  unloading  is  shown  in  Figure  4, 
where  the  bending  stress  by  elastic  deformation  is  illustrated  in  each  material. 


Figure  4.  Changes  in  stress  distribution  in  cross  section 
according  to  load  (Case  1) 


Case  2:  The  stress  in  a  ferromagnetic  layer  reaches  the  yield  stress,  after  SMA  layer 
reaching  the  transformation  stress  in  some  part. 

The  stress  distribution  of  Case  2  upon  loading  and  unloading  is  shown  in  Figure  5. 
Under  increasing  bending  moment  first  elastic  stress  distribution  (a),  then  the  stress  in  the 
SMA  layer  reaches  the  transformation  stress  Oo  at  the  position  of  yi  (b),  and  when  the 
transformation  domain  advances  to  yi=Yj,  a  ferromagnetic  layer  reaches  the  yield 
stress  Of  (c).  It  is  noted  in  (b)-(e)  that  Y]  remains  constant  until  y 3  reaches  Yj.  During 
unloading,  the  stress  decreases  first  in  elastically  in  all  domains  (d),  next,  the  stress 
becomes  constant  from  the  upper  part  of  the  SMA  layer  to  the  position  of  y3  where  the 
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stress  reached  reverse  transformation  stress  o i  (e).  In  addition,  after  the  stress  at  location 
y3=Yi  reaches  c h  the  stress  inside  portion  (y<y2)  decreases  elastically  (f).  Finally,  the 
stress  in  the  entire  SMA  layer  decreases  elastically  when  the  stress  in  the  SMA  on  the  top 
surface  becomes  small  than  Ci  (g). 


Unloading 

Figure  5.  Changes  in  stress  distribution  in  cross  section 
according  to  load  (Case  2) 


Case  3:  The  stress  in  a  ferromagnetic  layer  reaches  the  yield  stress,  after  the  entire 
domain  of  the  superelastic  SMA  layer  reaching  the  transformation  stress  Co. 

The  stress  distribution  of  Case  3  upon  the  loading  and  the  unloading  is  shown  in 
Figure  6.  In  early  stage  of  loading,  the  stress  in  a  ferromagnetic  layer  does  not  reach  the 
yield  stress  yet  even  after  the  stress  in  all  domains  of  the  SMA  layer  reaches  the 
transformation  stress  Co  (c).  A  neutral  axis  position  changes  with  an  increase  in  the  load, 
and  the  stress  reaches  the  yield  stress  of  finally  in  a  ferromagnetic  layer  (d).  The  process 
of  unloading  is  shown  in  Figure  6  (e)  to  (h). 
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Figure  6.  Changes  in  stress  distribution  in  cross  section 
according  to  load  (Case  3) 


For  each  stress  distribution  ox(y)  of  the  three  cases,  the  following  equations  are  valid, 
i.e.  the  equilibrium  of  force  and  moment. 


|V;r(y)3-rfy  =  0 

M  =-\hax{y)ybdy 


(1) 

(2) 


The  neutral  axis  position  and  the  relation  between  the  moment  and  the  curvature  are 
obtained  by  solving  these  equations.  Let  us  focus  on  the  case  2,  particularly  the  stress 
state  of  figure  5(b).  When  a  neutral  axis  position  is  and  the  transformation  stress 
position  is  yi,  the  stress  distribution  in  each  domain  becomes 
in  ferromagnetic  layer  (Ckychf) 

O-t-y)  =  £  JLzZ  (3) 

p 

in  SMA  layer  below  the  transformation  stress  Oo  (hf<y<yO 

o{y)  =  Esm^  (4) 

P 

in  the  transformation  domain  of  SAM  (yi<y<h) 

o(y)  =  -oa  (5) 
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By  substituting  Eqs.  (3),  (4),  and  (5)  to  Eqs.  (1)  and  (2),  unknown  ^  and  yi  are  solved 
and  they  are  given  by 


2l=4+_Sl£  (7) 

h  h  ESUAh 

Moreover,  by  substituting  Eqs.  (3)-(7)  to  Eq.  (2),  the  relation  between  the  normalized 
bending  moment  and  curvature  is  obtained  as 


Eq.  (8)  is  valid  for  the  range  of  curvature,  i.e.  from  the  curvature  with  transformation 
stress  Oo  in  top  :(y=h)  of  SMA  layer  to  the  curvature  with  yield  stress  Of  at  bottom  (y=0) 
of  ferromagnetic  layer.  This  range  of  the  curvature  is  given  by 


Similarly,  the  relations  between  the  bending  moment  and  the  curvature  for  the  three 
cases  of  Figs.  4-6  can  be  calculated.  The  results  for  cases  1,  2  and  3  are  shown  in 
APPENDIX.  The  conditions  under  which  three  cases  are  valid,  are  obtained  as 

Case  1 
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Case  2 


(12) 

(13) 


The  maximum  normalized  curvatures  in  these  cases  are  given  by 


Case  1  Case  2  Case  3 


The  maximum  deformability  of  the  composite  plate  can  be  analyzed  for  a  given  set  of 
the  mechanical  properties  and  the  thickness  ratio  of  materials  by  using  equation  (14). 


2.2.  Analytical  Results  and  Discussion  on  laminated  composite  for  bending 

The  relation  between  the  bending  moment  and  the  curvature  is  predicted  by  the 
present  model  for  two  type  of  the  composite,  i.e.  Fe/CuAIMn  and  FeCoV/CuAlMn. 
Figure  7(a)  is  the  idealized  stress-strain  curves  of  Fe  and  CuAlMn.  The  results  of  the 
predicted  relation  between  the  normalized  bending  moment  and  the  normalized  curvature 
for  thickness  ratio  hf/h=0.5  are  shown  in  Figure  7(b).  The  state  of  the  stress  for  this  case 
corresponds  to  Case  1,  Fig.4,  i.e.  the  stress  in  SMA  layer  is  not  superelastic  plateau,  thus, 
the  superelastic  loop  is  not  observed  as  evidenced  in  Figure  7(b).  Therefore,  the 
composite  plate  of  Fe  and  CuAlMn  is  undesirable  as  effective  bending  actuator 
component. 
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Figure  7.  (a)  Stress-strain  curve  for  Fe  and  (b)  Relation  between  normalized  bending  moment  and 
normalized  curvature  for  CuAlMn-Fe  composite. 


Next,  the  FeCoV/CuAIMn  composite  plate  was  analyzed  by  using  the  mechanical 
property  data  shown  in  Figure  8(a).  Figure  8(b)  shows  the  analytical  results  for  hf/h=0.5, 
exhibiting  clearly  superelastic  behavior.  By  using  FeCoV  whose  yield  stress  is  larger  than 
Fe,  yet  its  soft  magnetic  property  is  better  than  Fe,  we  can  achieve  now  the  state  where 
most  of  the  CuAIMn  layer  becomes  a  transformation  domain,  corresponding  to  almost 
the  state  of  Case  3.  Moreover,  the  maximum  curvature  was  2.22  times  larger  and  the 
bending  moment  was  1.60  times  larger  than  those  of  the  composite  with  Fe.  Therefore, 
the  FSMA  composite  so  identified  is  promising  as  an  effective  bending  actuator 
component. 


Figure  8.  (a)  Stress-strain  curve  for  FeCoV,  (b)  Relation  between  normalized  bending  mo  and 
normalized  curvature  for  CuAlMn-FeCoV  composite. 
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Next,  we  performed  a  set  of  parametric  studies  to  examine  the  effects  of  material 
parameters  (Of,  Ef,  Oo,  Oi,  Esma)  and  geometrical  parameter,  i.e.,  thickness  ratio  (hf/h). 
The  predicted  results  are  shown  in  Figure  8,  where  (a)-(f)  denote  the  case  of  changing 
parameters,  yield  stress  of  ferromagnetic  material  (Of),  the  upper  plateau  stress  (Oo)  and 
lower  plateau  stress  (Oi)  of  CuAIMn  superelastic  loop,  and  ratio  of  ferromagnetic  plate 
(hf)  to  the  composite  (h),  hf/h,  Young's  modulus  of  ferromagnetic  material  (Ef)  and  that  of 
SMA  (Esma),  respectively. 


(a)  a,  (b)  o0 


Figure  9.  Change  in  superelastic  behavior  of  bending  plate  influenced  by  various  parameters, 

(a)  yield  stress  of  Fe,  (b)  upper  transformation  stress  of  SMA,  (c)  lower  transformation  stress  of  SMA,  (d) 
thickness  ration  of  Fe  to  FSMA  composite,  (e)  Young's  modulus  of  Fe  and  (f)  Young's  modulus  of  SMA. 

When  the  yield  stress  of  the  ferromagnetic  material  increases,  it  is  clear  from  Figure 
9(a)  that  both  bending  moment  and  the  curvature  increase.  When  transformation  stress  Co 
of  SMA  increases,  it  is  found  from  Figure  9(b)  that  the  bending  moment  increases  and  the 
curvature  decreases.  It  can  be  seen  from  Figure  9(c),  the  lower  limit  of  the  superelastic 
loop  decreases  if  the  reverse  transformation  stress  Oi  decreases.  When  the  thickness  of 
the  ferromagnetic  layer  increases,  it  is  clear  from  Figure  9(d)  that  the  bending  moment 
increases  though  the  curvature  decreases.  Oppositely,  because  the  thickness  of 
superelastic  SMA  layer  increases  when  the  thickness  of  a  ferromagnetic  layer  decreases, 
the  superelasticity  behavior  increases.  Therefore,  the  bending  moment  decreases,  and  the 
curvature  increases.  From  Figure  9(e),  the  maximum  curvature  decreases  though  the 
bending  moment  does  not  change  when  the  Young's  modulus  of  the  ferromagnetic 
material  increases.  Therefore,  an  increase  in  the  Young's  modulus  of  the  ferromagnetic 
material  is  undesirable  as  the  composite.  From  Figure  9(f),  the  bending  moment 
decreases  when  the  Young's  modulus  of  SMA  increases.  The  design  of  a  more  high 
performance  ferromagnetic  shape  memory  alloy  composites  becomes  possible  by  the 
materials  design  based  on  the  above  analysis. 


13 


2.3.  Modeling  on  superelastic  behavior  of  coil  spring  made  of  a  composite  wire  with 
rectangular  cross  section 

With  the  aim  of  designing  a  high-speed  linear  actuator,  the  superelastic  characteristic 
of  a  coiled  spring  of  the  ferromagnetic  shape  memory  composite  wire  with  rectangular 
section  is  analyzed.  Figure  9  shows  the  analytical  model.  The  magnetic  force  is  generated 
in  the  ferromagnetic  material  by  the  magnetic  field  gradient,  and  displacement  is 
generated  in  the  spring  by  the  hybrid  mechanism  described  in  Introduction.  The  relation 
between  this  spring  force  and  displacement  is  analyzed. 


D:  the  diameter  of  spring  (D=2R) 

d:  the  diameter  of  wire 

p:  the  pitch  of  one  cycle 

n:  the  number  of  turns 

L:  the  length  of  spring  without  load  (L=np) 

a:  the  inclined  angle  of  the  wire  to  the  x-y  plane 

Figure  10.  Analytical  model  of  coil  spring  with  rectangular  cross  section. 

When  axial  force  P  is  given  to  the  spring,  the  wire  of  the  ferromagnetic  shape 
memory  composite  is  subjected  to  torque  T.  The  relation  between  spring  force  P  and 
torque  T  is  given  by  the  next  equation  by 

T  =  PRcosa  (15) 
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for  a  twist  angle  per  unit  length  of  the  rectangular  section  wire  of  CD,  the  total  twist  angle 
<j>  is  2njiRcoseca  as  the  total  length  of  the  wire  is  2n7iRseca.  Therefore,  the  displacement 
of  the  spring  is  calculated  by  the  next  equation. 

S=S  .  +8 .  sS 

torsion  shear  torsion 

=  R<p  =  2nnR^meca  (16) 

It  is  assumed  in  the  present  model  that  the  displacement  due  to  direct  shear,  Sshear  is 
neglected.  This  is  justified  for  large  ratio  of  D  to  a  or  b.  Then,  the  relation  between  the 
,  spring  force,  P  and  displacement,  8  can  be  calculated  if  the  relation  between  the  twist 
angle  per  unit  length  (0  and  the  torque  T  of  the  rectangular  section  wire  is  known,  which 
will  be  obtained  in  the  following. 

Analytical  Model  for  Torsion  of  Composite  Wire  with  Rectangular  Section 

To  generate  large  magnetic  force  by  the  hybrid  mechanism,  it  is  necessary  to  increase 
the  area  of  a  ferromagnetic  material  in  the  rectangular  section,  while  meeting  the 
requirement  that  the  ferromagnetic  material  should  not  reach  its  yield  stress.  The  stress 
field  in  the  rectangular  section  can  be  calculated  from  the  shear  strain  distribution  of  the 
rectangular  section  for  a  given  twist  angle. 

Let  us  look  at  the  rectangular  section  of  a  composite  with  width  2a  and  height  2b  as 
showing  in  Figure  11.  We  introduce  the  assumption  that  the  spring  deformation  is 
uniform  along  the  wire  direction  (z-axis)  and  plane  displacements  u  and  v  are  in 
proportion  to  z,  as  follows; 

u=-oyz,  v=  cacz,  w -oxp^x,y)  (17) 

where  the  function  <p(x,y)  is  the  Saint- Venant's  function  [11]  that  satisfies  the  equilibrium 
equation  and  2D  compatibility  equation  of  strain.  For  the  spring  with  rectangular  cross 
section,  the  shear  strain  components  are  expressed  as 
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(c)  Y/aco 


Figure  1 1 .  Contour  line  distributions  of  shear  strain  in  rectangular  section, 

(a)  Yzx/aco  ,(b)  Yz/ao)  and  (c)  y/m,  where  a  is  the  length  of  longer  side  of  a  rectangular  cross  section  of 
a  FSMA  composite  and  co  is  the  twist  angle  per  unit  length. 


(b)  Yz/aco 
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(18) 

(19) 


Therefore,  the  effective  shear  strain  acting  on  the  rectangular  cross  section,  yis  calculated 

by 


r^yjrJ+rJ  (2°) 


For  a=2  and  b=l,  the  contour  line  distributions  of  shear  strain  component  yzx,  yzy  and 
effective  shear  strain  y  divided  by  act)  are  shown  in  Figure  11(a),  (b)  and  (c)  respectively. 
y7X  becomes  0  at  x=-a  and  a,  and  it  reaches  to  the  minimum  value  at  y=b  on  the  y  axis, 
and  becomes  the  maximum  at  y=-b  on  the  y  axis.  yzy  reaches  to  the  minimum  value  at  x=- 
2,  y=0,  and  becomes  the  maximum  at  x=2,  y=0.  The  normalized  effective  shear  strain, 
y/aco  reaches  the  maximum  value  0.930  at  the  center  of  long  side  edges,  and  reduces 
toward  the  center. 

The  effective  shear  stress  induced  in  the  ferromagnetic  material  is  calculated  by 
multiplying  yby  the  shear  modulus  Gf  of  the  ferromagnetic  material.  The  effective  shear 
stress  distribution  of  the  ferromagnetic  material  in  the  rectangular  section  is  calculated  for 
a  given  set  of  twist  angle  per  unit  length  (0,  size  a  and  b.  Then,  the  optimum  shape  of  the 
ferromagnetic  material  can  be  determined  from  its  domain  under  the  condition  that  the 
effective  shear  stress  does  not  exceed  the  yield  stress  in  shear  tf  of  the  ferromagnetic 
material. 

If  FeCoY  (Gf=70GPa,  Xf=231MPa)  is  used  as  a  ferromagnetic  material,  and  CuAIMn 
is  used  as  a  superelastic  SMA,  then  for  co=0.003,  a=2,  and  b=l,  y/{oa<0.55  is  obtained 
from  the  requirement  of  Gf  y<tf.  Figure  1 1  shows  the  optimized  rectangular  section  of  the 
composite  obtained  by  this  design,  where  the  dark  area  of  FeCoV  satisfies  y/coa  <  0.55. 
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Figure  12.  Optimized  rectangular  section  of  the  composite  made  of  ferromagnetic,  FeCoV  and 

superelastic  SMA,  CuAlMn. 

Next,  we  examine  the  relation  between  the  twist  angle  per  unit  length  to  and  the 
torque  T  of  the  composite  wire  with  rectangular  section.  The  torque  is  calculated  by 

r  =  Mz=/*«cosa=J  M  (21) 

=l 

Here,  we  can  define  three  domains  in  the  composite  during  loading. 

Domain  1:  Domain  of  ferromagnetic  material 

Domain  2:  Domain  with  effective  shear  stress  less  than  the  forward  transformation 
shear  stress  of  SMA,  To 
Domain  3:  Transformation  domain  of  SMA 
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The  effective  stress  in  the  ferromagnetic  material  is  obtained  by  multiplying  shear 
modulus  Gf  by  the  corresponding  effective  shear  strain  for  the  elastic  deformation.  In 
SMA,  it  is  necessary  to  judge  if  the  effective  shear  stress  X  is  below  the  forward 
transformation  shear  stress  %  x  is  obtained  by  multiplying  shear  modulus  Gsma  by  y  if 
the  effective  shear  stress  of  domain  2  is  below  the  forward  transformation  shear  stress  To- 
In  domain  3  where  the  effective  shear  stress  X  reaches  the  upper  transformation  shear 
stress  to,  then  x  =  to- > 

Because  a  shear  strain  component  proportionally  increases  with  an  increase  in  to,  by 
multiplying  the  corresponding  shear  strain  component  by  the  modified  shear  modulus 
G=to/y,  the  shear  stress  component  for  which  x  becomes  to  is  calculated.  That  is. 
Equation  (21)  is  applicable  to  domain  3  by  using  Equation  (24).  Then,  torque  T 
corresponding  to  the  twist  angle  per  unit  length  to  is  calculated  by  equation  (21)  by  using 
the  modified  shear  modulus  in  each  domain  according  to  the  following  Eqs. 


(22) 


(23) 


(24) 

where,  (Of  is  input  data,  and  it  is  the  maximum  twist  angle  per  unit  length  when  the  cross 
section  is  optimized,  (Of  =0.0015  for  cross  section  shape  of  Figure  12. 


Domain  1:  ...  T/(0 

G/eof 


G  =  G, 


Domain  2:  ^  T/w  and  ,  ro 

Gf®J  Gsm 


G  =  G,. 


Domain  3:  T/w  and  Y>Jh_ 

G/®/  Gsma 


G 

y 
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Next,  the  case  of  unloading  is  considered.  The  stress  in  each  domain  decreases  during 
the  unloading,  but  the  superelasticity  in  SMA  was  generated  in  domain  3  where  the 
effective  shear  stress  reached  the  transformation  stress  during  the  preceding  loading,  it  is 
necessary  to  divide  domain  3  into  three  sub-domains. 

Domain  3-1:  x  above  the  reverse  transformation  stress  Xi 
Domain  3-2:  x  equal  to  the  reverse  transformation  stress  Xi 
Domain  3-3:  x  beldw  the  reverse  transformation  stress  X] 

For  domain  3-1,  the  effective  shear  stress  x  is  larger  than  the  reverse  transformation 
stress  Xj.  The  shear  stress  component  of  x=Xo  is  calculated  by  multiplying  modified  shear 
modulus  of  Equation  (24)  by  the  shear  strain  component,  and  it  decreases  from  this  stress 
state  elastically  in  proportion  to  Gsma  in  domain  3.  That  is,  the  shear  stress  component  is 
calculated  by  multiplying  the  modified  shear  modulus  of  Equation  (26)  by  the  shear 
strain  component  in  the  range  of  the  effective  shear  strain  of  Equation  (25). 

Domain  3-1  y>_JL_  and  fa-r,>a  (25) 

Gjam  (o}{  -  w) 

G^~Gstu  (26) 

For  domain  3-2,  because  the  effective  shear  stress  reaches  the  reverse  transformation 
stress  Xi,  the  shear  stress  remains  constant,  i.e.  x=Xi.  That  is,  the  shear  stress  is  calculated 
by  multiplying  the  modified  shear  modulus  of  Equation  (28)  by  the  shear  strain  in  the 
range  of  effective  shear  strain  of  Equation  (27). 

Domain  3-2  ^  and  (27) 

Gsma  Gsua  (ty  —  (o) 

G=—  (28) 

7 
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For  domain  3-3,  the  superelasticity  disappears  because  the  effective  shear  stress 
lowers  more  than  The  range  of  effective  shear  strain  and  modified  shear  modulus  are 
given  by 

Domain  3-3  (29) 

G  =  Gsma  (30) 

i 

The  torque  T  corresponding  to  to  can  be  analyzed  from  equation  (21)  by  calculating 
effective  shear  strain  of  each  area  using  the  modified  shear  modulus  corresponding  to 
each  domain  defined  by  Equations  (22),  (23),  (24),  (26),  (28)  and  (30).  The  relation 
between  the  force  and  displacement  of  a  spring  can  be  calculated  by  using  Equation  (15) 
and  (16). 

2.4.  Analytical  results  and  discussion  on  composite  spring 

Based  on  the  above  model,  we  made  predictions  of  the  torque  (T)  -  twist  angle  (to) 
relation,  and  also  of  the  spring  force  (P)  -  displacement  (5)  relation  where  the  idealized 
stress-strain  relations  of  ferromagnetic  FeCoV  and  superelastic  CuAIMn  shown  in  Figure 
13  are  used. 


Shearing  strain :  y 

Figure  13.  Idealized  stress-strain  curves  of  FeCoV  and  CuAIMn. 
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(a)  Relation  between  torque  and  (b)  Spring  force  (P)-displacement  (8)  curve 

normalized  twist  angle 

Figure  14.  Superelastic  behavior  of  Fe/CuAIMn  composite  spring. 

Figure  14  show  the  analytical  results  for  the  case  of  maximum  twist  angle  per  unit 
length  w=0.003  of  a  composite  plate  wire  with  a=2mm  (width  is  4mm),  and  b=lmm 
(height  is  2mm).  Figure  14(a)  shows  the  relation  between  the  torque  and  the  normalized 
twist  angle,  indicating  that  the  torque  rises  proportionally  as  the  twist  angle  increases,  and 
the  transformation  of  SMA  begins  at  O>a=0.0025,  reaching  the  transformation  stress  with 
coa=0.0042  in  all  domain  of  SMA.  After  (Da  reaches  0.006,  the  superelastic  loop  exhibits 
the  reverse  transformation  corresponding  to  the  unloading. 

Figure  14(b)  shows  the  relation  between  the  spring  force  and  the  displacement  of  the 
coil  spring  of  length  L=  100mm,  diameter  D=25mm,  pitch  p=5mm  and  number  of  turns 
n=20.  The  maximum  displacement  of  this  coiled  spring  was  59.2mm,  the  spring  force 
became  78.4N. 

We  made  a  parametric  study  to  examine  the  effects  of  each  parameter  on  the  P-8 
relation.  Figure  15  shows  the  analytical  results  of  the  P-5  relations  influenced  by  various 
parameters,  (a)  Gsma.  (b)  To  ,  (c)  Gf,  (d)  if  and  (e)  Ti.  From  Figure  15  (a),  it  is  clear  that 
shear  modulus  of  superelasticity  SMA  does  not  influence  the  maximum  displacement  and 
the  maximum  spring  force.  It  is  noted  from  Figure  15  (b),  that  the  spring  force  increases 
with  an  increase  in  forward  transformation  shear  stress  %  It  is  clear  from  Figure  15  (c), 
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that  the  spring  force  does  not  change  and  only  the  maximum  displacement  increases  if  the 
shear  modulus  of  the  ferromagnetic  material  becomes  small  resulting  in  larger 
displacement  of  the  spring.  It  can  be  seen  from  Figure  15  (d),  that  both  the  spring  force 
and  displacement  increase  the  superelastic  behavior  when  the  yield  stress  of  the 
ferromagnetic  material  increases.  It  is  noted  from  Figure  15  (e),  that  the  lower  limit  of 
superelastic  loop  decreases  if  the  reverse  transformation  stress  Xi  decreases. 

I 

In  summary,  lager  Xf  of  the  ferromagnetic  material  and  softer  ferromagnetic  material 
will  provide  a  spring  actuator  with  larger  displacement.  And,  to  obtain  large  force  of  the 


spring,  use  of  SMA  of  larger  Xo  is  desired. 


(a)  Gsma 


(c)Gf 


(b)To 


(d)tf 
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100 


(e)ti 

Figure  15.  Effects  of  various  parameters  on  P-8  relation  of  FMSA  composite  springs, 
(a)  SMA  shear  modulus,  Gsma  >(b)  forward  transformation  shear  stress,  To, 

(c)  shear  modulus  of  a  ferromagnetic  material,  Gf,  (d)  the  yield  stress  in  shear  of  a 
ferromagnetic  materia!, Tf,  and  (e)  reverse  transformation  shear  stress,  Tj. 


(a)  Shape  of  cross  section 


(b)  Spring  force-displacement  curve 


Figure  16  Superelastic  behavior  of  Fe/CuAIMn  composites. 


Next,  we  shall  compare  the  mechanical  performance  (P-8  relation)  of  a  spring 
between  "rectangular"  and  "square"  cross  section.  To  this  end,  the  cross  section  area  of 
the  square  is  made  equal  to  that  of  the  rectangular  studied  earlier  (Figure  13).  The 
analytical  results  of  the  optimum  square  cross  section  of  FeCoV/CuAIMn  composite  are 
shown  in  Figure  16(a),  while  the  P-5  relation  of  the  FSMA  spring  with  this  square  cross 
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section  is  given  in  Figure  16(b)  as  a  dashed  line  where  the  results  of  the  rectangular  cross 
section  are  also  shown  by  solid  line.  A  comparison  between  the  square  cross  section  of 
Figure  16(a)  and  the  rectangular  cross  section  of  Figure  12  reveals  that  the  FSMA 
composite  spring  with  square  cross  section  provides  larger  force  capability  than  that  with 
the  rectangular  cross  section  for  the  same  cross  section  area.  However,  the  effectiveness 
of  using  the  spring  with  the  square  cross  section  remains  to  be  determined  after  its 
effectiveness  of  inducing  large  magnetic  force  between  the  neighboring  turns  of  the 
spring. 

2.5.  Concluding  remarks 

The  predicted  results  of  the  bending  moment  -  curvature  of  a  FSMA  composite  plate 
exhibit  superelastic  behavior  of  the  composite  beam  while  those  of  the  FSMA  composite 
spring  with  rectangular  cross  section  show  also  similar  superelastic  behavior.  The  above 
superelastic  behavior  is  the  performance  required  for  FSMA  composite  actuators  with 
high  force  and  displacement  capability.  The  results  of  the  simple  model  were  used 
effectively  for  optimization  of  the  cross  section  geometry  of  two  types  of  FSMA 
composite,  bending  and  torsion  types. 
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Chapter  3.  Modeling  of  Effective  Magnetic  and  Mechanical  Properties  of 

Fe-NiTi  Particulate  Composite 

In  this  chapter,  we  will  discuss  two  models  to  predict  the  magnetic  properties  of 
FSMA  particulate  composites,  first  for  magnetic  permeability  and  second  for  saturation 
magnetization. 

3.1.  Composite  Magnetic  Permeability 


x,  x, 

Actual  composite  Eshelby’s  model 


Figure  17.  Eshelby’s  model  for  a  FSMA  particulate  composite 
In  the  absence  of  inhomogeneities,  see  Figure  17, 

B  =  (31) 

Here  B  is  the  magnetic  flux,  Hm  is  the  magnetic  field  intensity  in  the  matrix  and  jum  is 
the  magnetic  permeability  of  the  matrix.  As  shown  in  figure  17,  the  matrix  phase  (m- 
phase)  is  the  non-ferrous  NiTi  alloy  (paramagnetic  material)  and  Iron  in  particulate  form 
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constitutes  the  fiber  phase  (f-phase).  Since  the  magnetic  susceptibility  %  for  paramagnetic 
materials  typically  lies  in  the  range  of  10"3-  10'5,  the  magnetization  vector  is  very  weak. 

Therefore,  equation  (31)  takes  the  following  form  [12]:- 

Hm  =  Ha  +M  <t=>  Hm  «  Ha 

M  ~  0  — »  Magnetization  vector  for  paramagnetics  (NiTi)  (32) 

B  =  Ho  . 

This  error  is  found  to  be  of  the  order  of  0.01%  which  is  considered  acceptable  by  the 
PI.  Upon  addition  of  inhomogeneity  (ferromagnetic  phase),  the  flux  and  field  vectors  are 
related  as  follows  in  Q;- 

B  +  B,,  =fi/'(H0  +HJ)  =  fim -(H0  +Hd-H* )  (33) 

By  definition  for  ferromagnetics, 

B  =  fj0-(H0+M-N-M)  =  fir(H-N-M)  (34) 

where  N  — >  demagnetization  factor,  a  tensor  quantity  depending  only  on  the  shape  of  the 
inhomogeneity,  analogous  to  Eshelby’s  tensor  of  mechanical  case. 

Choosing  Hd  =  ~N  ■  M  in  equation  (34)  we  get, 

Hd  =N  H*  =Sm  H* 

°r, 

B  +  Bd=nr(H0+Sm-H')  =  vm-(H0+Sm -H'-H')  { ^ 

Sm  =  N  "Eshelby's  Tensor"  for  Magnetic  case 
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From  equation  (35)  H*  can  be  solved.  When  the  volume  fraction  of  ferromagnetic 
inhomogeneties  if)  becomes  finite,  the  interactions  between  ferromagnetic  particles  need 
to  be  accounted  for.  Mori-Tanaka’s  mean-field  theory  has  been  used  in  the  present  case1. 
We  thus  define  H  as  the  average  magnetic  field  and  equation  (35)  gets  modified  to:- 

B  +  Bd=nf-(H0+H+Hd)  =  vm-(H0  +  H+Hd-H') 

or,  (36) 

Since  the  volume  integration  of  Hd  must  vanish  over  the  entire  domain,  we  get, 
H=-f(Hd  -//’)  (37) 

Substituting  in  equation  (36)  and  simplifying  gives  us, 

Following  the  steps  enumerated  in  section  4.2.1  of  [13]  the  concentration  factor 
tensor  A  is  evaluated, 

Hf  =  Ha  +  H  +  Hd  =  H„  +  Hd(\ - f)  +  fH'  =  Ha  +  \  1  -  j)Sm  +  fl\ //* 

or,  (39) 

Or, 


1  Atomic  interaction  effects  have  been  neglected  here.  This  is  something  that  the  PI  will  be  considering  in 
future. 

2  Section  on  “Effective  Property  tensor  of  a  composite:  Field  vector  (Z0 )  boundary  condition” 
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(40) 


or, 

Hf=A-H0 

with  A  =  [l  +  il-f)Sm+fl}-hf-/Jm)fa-f)Sm+fl}+  M.Y  -fa-ltf)] 

Thus, 

f*c  =Mm+f(Vf  -J*r*y-A  (41) 

juc  —»  Composite  magnetic  permeability. 

Thus,  by  knowing  the  properties  of  the  individual  constituents  of  the  composite  and  the 
geometry  of  the  inclusion,  the  composite  magnetic  permeability  can  be  evaluated. 

3.2.  Saturation  Magnetization  (Mcs)  of  Fe-NiTi  particulate  composite 

The  magnetic  field  intensity  in  the  composite  is  evaluated  as  follows3:- 
Hc=(l-f)(Hm)+f(Hf) 

(Hf)  =  H0  +  H  +  Hd 

=  H0-f(Hd-H*)+Hd 
=  H0+(l~f)Hd+fH' 

(Hm)  =  H0  +  H 

=  H0-fHd+fH* 

hc  =  h0+m; 

=  HD+fH * 

=  H0-fMs  (42) 

or, 

Mcs=fMs 

Here,  Ms  ->  Saturation  Magnetization  of  the  ferro  -  particulate  (f  -  phase) 
and  Mcs  — >  Composite  Saturation  Magnetization. 


3  ^  Denotes  volume  average  over  the  entire  domain. 
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Since  the  experimental  data  tabulates  the  Saturation  magnetization  for  various  Fe 
concentrations,  the  composite  saturation  magnetization  is  a  very  useful  comparison 
parameter. 

We  measured  the  saturation  magnetization(Ms)  of  Fe-TiNi  particulate  composites 
that  we  processed  by  using  Spark  Plasma  Sintering(SPS),  see  chapter  4  where  four 
different  weight  percents  of  Fe  phase  are  used  to  form  Fe-TiNi  composites.  We  compare 
the  predicted  and  measured  Ms  values  in  Table  1,  and  also  in  Figure  18. 


Weight  %  Fe  (w) 

imiifiin 

Mcs  (Experimental) 

M[  (Predicted) 

30 

0.26133 

40 

52.266 

50 

0.4522 

93 

70 

0.6582 

135 

131.64 

100 

1.0 

200 

200 

Table  1.  Comparison  of  saturation  magnetization,  Ms  between  the  predictions  and  experimental  data 


Saturation  Magnetization  vs.  Volume  Fraction 


Figure  18.  Comparison  of  saturation  magnetization  predicted  by  our  model  and  measured 


4  /  = - ;  pFe ,  pNm  — >  densities  of  Fe  and  NiTi  (7874  and  6500  g/cc  respectively) 

1  I  1-W  pFe 

W  P  NiTi 
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It  is  clear  by  comparing  the  model  predictions  and  the  experimental  results  that  there 
is  not  a  100  %  match  in  the  values.  Especially,  for  lower  volume  fractions  the  model 
overestimates  the  Saturation  Magnetization  of  the  composite.  But  the  margin  of  error 
seems  acceptable.  The  model  could  be  further  improved  by  taking  into  account  the 
interaction  effects  which  the  PI  is  presently  investigating. 

3.3.  Stiffness  properties  of  a  FSMA  composite 

Here  we  constricted  an  analytical  model  to  predict  the  effective  stiffness  and  slope  of 
the  stress-strain  (SS)  curve  of  a  FSMA  particulate  composite  which  is  composed  of  the 
particulate  SMA  phase  and  ferromagnetic  matrix  phase.  The  model  is  restricted  to 
particulate  geometry,  and  also  the  role  of  particulate  and  matrix  can  be  exchanged.  The 
model  is  based  on  Eshelby’s  model.  A  typical  SS  curve  of  a  composite  containing  SMA 
phase  is  idealized  as  in  Fig.  19  which  consists  of  three  linearized  stages,  the  first  stage 
(100%  austenite  phase,  E*m ),  second  stage  a  mixture  of  austenite  and  martensite  phases 
(ENiTi )  and  the  third  stage  (100%  martensite  phase,  E^1T‘ ).  Please  note  that  in  the  second 

stage,  the  fraction  of  martensite  phase  increases  with  increase  in  stress,  thus  this  stage  is 
often  called  “stress-induced-martensite  phase  transformation  (SIM)”.  SIM  is  a  key  part  of 
superelastic  behavior,  shown  as  the  loop  in  Fig.  19.  In  the  following,  we  state  in  details 
the  sequence  of  our  model  calculations. 


Figure  19.  Linearized  version  of  SS  Curve  of  a  SMA  composite  where  the  first  and  the  third 
stages  are  100%  austenite  and  100%  martensite  phase  respectively,  and  the  second  stage  is  a 
mixture  of  both  phases,  i.e.  stress-induced  martensite  phase  transformation  (superelastic  stage) 
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For  the  first  stage  of  the  loading  curve,  the  transformation  strain  £x  =  0 .  The  problem 
thus  reduces  to  Eshelby’s  inhomogeneity  problem  (Figure  20)  where  the  applied  stress 
causes  a  relative  strain  among  the  matrix  and  fiber.  Using  Mori-Tanaka’s  Mean  Field 
Theory  in  association  with  Eshelby’s  the  stresses  in  the  fiber  Q  domain  are  given  by:- 

O°+0  =  ENm-(£O  +  £  +  £)=EFe-(£o+£  +  £-£*)  (43) 

Here, 

o°  — >  applied  stress  field 
a  ->  disturbance  stress  field 
£->  disturbance  strain  field 
£  — >  average  strain  due  to  interactions 
EFe  ->  stiffness  tensor  of  Fe 


Cl  ->  Fiber  domain 
D-Q  —» Matrix  domain 


(a)  -(b) 


Figure  20.  Eshelby’s  Inclusion  problem  for  use  in  the  simulation  of  first  and  third  stages  behavior 

in  SS  curve  of  a  SMA  composite 
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By  definition, 

£#=(EFer‘o°  (44) 

or, 

<J  =  EFe  •(£  +  £-£*)  (45) 

The  volume  average  strain  disturbance  over  the  entire  domain  is  zero,  thus  we  get:- 
£  =  -/(£-£*)  (46) 

Substituting  equation  (46)  in  (43)  and  after  some  algebra  we  get:- 

Ec  -Vr+/Kf  (47) 

where 

£*=K  0° 

K  =  [(ENiT1  -EF')-((l-/)S  +  /l)+EFe]"1  -(ef*  -ENm)-(EF*r  (48) 

and  where 


ENiTi  ->  E  Am ,  during  stage  1  loading  and  ENiTl  -»  E^ITl ,  during  stage  3  loading 
Calculation  of  o„s 


From  equation  (47)  we  can  find  the  slope  of  the  compression  SS  curve  of  the 
composite  given  the  input  parameters,  namely,  material  properties  of  constituent  phases 
NiTi  and  Fe.  But  in  order  to  trace  the  SS  curve  of  the  composite,  we  need  to  evaluate  the 
stress  required  to  be  applied  to  the  specimen  such  that  the  NiTi  phase  of  the  composite 
has  a  stress  value  equal  to  Martensite  Start  Stress  (  o*1 ).  The  following  set  of  equations 


describes  a  way  of  doing  that:- 

°Ms  =  o°  =  V0<T33  =  (1  -  /)«Fe  +  /«m™  (49) 

oFe  =(I-/EFe  •  (S  - 1)  •  K  A  )  •  o°  =>  aFe  =  RA  •  o°  (50) 

where 

Vo  =  [001000f  ;RA=(I-EFe  /(S-I)  KA);  (51) 

K  A  =  [(EAiTi  -  EFe )•  {(1  -  /)S  +  /!}+  EFe  Y  ■  (EFe  -  EAiTi  )•  (EFe )_I  (52) 
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=*  (I  •  RA  )V„ff„  -  jsC ;  RAB  =  (I  -  RA  )v  =>  a„  =  (53) 

Thus,  applying  er33  stress  in  3-direction  would  produce  SIM  in  the  NiTi  phase  of  the 
FSMA  composite.  Similarly,  the  total  strain  in  the  composite  is  given  by:- 

(®)i total  =.(1-/X*)fc  +  /<e)Nm  ^(£)to,a.  =a-/)(£0+£)  +  /(S°+i  +  e) 

=»<«>« «•+/»’  =•<■=>,«,, =[VT +/*«}«• -[^-r+/K*]-v^„  <54) 

Thus,  the  point  (  <r33 ,  (c)  (3,1) )  indicates  the  first  “kink”  point  of  the  assumed 

1  linear  SS  curve  of  the  FSMA  particulate  composite. 


In  the  second  stage  of  loading,  an  eigenstrain  et  ,  phase  transformation  strain  from 
Austenite  to  Martensite,  exists  in  the  NiTi  fiber  region.  Since  the  eigen  strain  or  the  non¬ 
elastic  strain  exists  only  in  the  fiber  phase,  the  problem  is  identical  to  Eshelby’s 
inhomogeneity  inclusion  problem  (Figure  21).  This  can  be  represented  as  follows:- 


In  Q,  (s°  +  e  +  e-et)=Ef*  -(e0  +e  +  e-£T  -£*)=EFe  (e"  +E  +  E-E**)  (55) 

where 

£**  =e*+£T,et  =[~vNm  ~vNm  1  0  0  Of -4 

£**  =  V  •  o°  +  W  •  et 

V  =  tei-EFe)-{(l-/)S  +  /l}+EFel'^(EFe-E^O-(EFer  (56) 

W  =  te‘-EF*>{(1-/)S  +  /I}+EF*J-1-E^i 


Since  the  SS  curve  of  pure  NiTi  is  assumed  to  be  linear,  E^m ,  the  stiffness  tensor  of  NiTi 
phase  exhibiting  et  eigen  strain  is  taken  as 


where  sj0tal  is  the  total  transformation  strain  during  compression  loading  of  NiTi. 
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The  total  volume  averaged  stress  in  domain  D  becomes 

(°)D  =  o°  =  f  1  -  /)EFe  •  (£°  +£)  + f  (a) { 

=>  «#  =  ( 1  -  •  (S°  - /( s  -I) •£**)  +  /{o)f 

=>  o°  =  ( i  -  /)<**  -  fa  -  mTe  •  <s  - 1)  •  s* + /{o)f 

=>  [i + ( 1  -  y)EFe  •  (S  - 1)  ■  ■ v]'  o°  +  ( 1  -  f) EFe  ■  (S  - 1)  ■  W  •  £t  =  (o)f 


[i+ri-/)EF'(s-Dv]- 


ff33+(l-/)Ete-(S-I)-W 


^NiTi 
—  V  NiTi 

1 

0 

0 

0 


4=K  (58) 


/?(6,1)(T33  +  7(6, 1)4  =  (o)r 


o°  x3  c° 


(a)  (b) 

Figure  7.  Eshelby’s  Inhomogeneity  Inclusion  problem  for  use  in  the  simulation  of  the  second 
stage  behavior  of  SS  curve  of  a  SMA  composite. 
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^  i?(3,l)<733  +  T(3,l)£33  —  (c33)  j  C33 
where 


(g33)/-mi)4 

/?(3,1) 


'O' 

VNiTi 

0 

~  V  NiTi 

R  =  [l+(l-/)EFe-(S-I)-v]- 

1 

0 

and  T  =  (1  -  /)EFe  •  (S  - 1)  •  W  • 

1 

0 

0 

0 

0 

0 

(59) 


The  volume  averaged  strain  is  computed  as:- 

(t)D=(l- f)(t°  +s)+  f(t°  +£+£-et)  =  z°  +  f{t~  -£t) 

=>  (e)d  =  |ef'  y  +  yv)-  a0  +  /( W  - 1)  ■  sT  i 

where 


'o' 

VNiTi 

0 

~VNiTi 

1 

0 

<733,£T  = 

1 

0 

0 

■ 

0 

0 

0 

(60) 


(61) 


Thus,  the  applied  stress  cr33  required  to  cause  any  given  e33  transformation  strain  can 
be  calculated  using  equation  (57).  Consequently,  the  total  strain  exhibited  by  the 
composite  for  the  applied  stress  can  be  evaluated  using  (58).  Plotting  cr33  versus  (e)d(3, 

1)  values  obtained  by  varying  e%3 between  (0,ejolal)  gives  us  the  SS  curve  for  the  Fe-NiTi. 
It  is  to  be  noted  that  equation  (57)  makes  use  of  ‘older’  e33to  compute E^1  (4) •  An 
iterative  process  for  the  calculation  of  cr33  has  been  proposed  here.  e33  is  varied  between 
(0,  ej0tal )  and  at  each  nth  step,  E^m(£33)  is  calculated  using  By  increasing  the 

number  of  steps,  the  error  can  be  minimized.  The  computations  are  performed  using 
MATLAB. 


36 


c)  3rd  stase  loading  curve 


Third  stage  loading  curve  is  similar  to  the  first  stage  in  the  sense  that  there  does  not 
exist  any  transformation  strain  in  the  NiTi  phase.  Thus,  the  problem  can  be  solved  as 
Eshelby’s  inhomogeneity  (Figure  20).  The  first  step  involves  calculation  of  martensite 


finish  stress  o£,f  as  enumerated  overleaf:- 


°Mf=  °0  =  ff33Vo  =(l  -/Ve  +  /<C 

oFe  =  (I  -EFe  •  /(S  -I)  •  Km ) •  o°  =>  cFe  =  RM  •  o° 
where 

Vo  =  [001000r;RM  =  (I-EFe-/(S-I)KM); 

Km  =  [(E»m  - EFe  )•  {(1  -  /) S  +  fl}+  EFe  Y  •  (e*  - ES™  )•  (EFe  Y 
=>  (l  *  RM)Vo£r33  =  /<C‘ ;  RMB  =  (I  -  RM)V0  =>  <r33  = 


(62) 

(63) 

(64) 

(65) 

(66) 


Thus,  applying  a33  stress  in  3-direction  would  produce  o„^(3,l)  in  NiTi  phase  of  the 
FSMA  composite. 

The  total  strain  in  the  composite  is  given  by:- 


=V-  fi£) F«  +  /(£)NiTi  =>  (8Ll  =(J- /X£°  +  £)  +  +C  +  f-) 

=><*), *.,  =  £°  +  ^  =><«)„„.  =  l(EFe)_1  +yKM]-o°  =  [(EF«r  +  /KM].V0a; 


(67) 
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Thus,  the  point  (  cr33>(£),ora;  (3,1))  indicates  the  second  “kink”  point  of  the  assumed 
linear  SS  curve  of  the  FSMA  particulate  composite. 
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Chapter  4.  Processing  of  Particulate  FSMA  Composites 

One  of  the  proposed  FSMA  composite  in  the  present  work  is  a  “particulate 
composite”  composed  of  particulate  powders  of  SMA  (NiTi)  of  superelastic  (SE)  grade 
and  ferromagnetic  material,  i.e.,  soft  iron  (Fe).  The  ordinary  metallurgical  route  for 
processing  particulate  composites  using  powders  i.e.,  standard  sintering  with  long 
duration  of  high  temperature  produces  unwanted  reaction  products  destroying  the 
original  properties  of  SMA  and  ferromagnetic  material.  Therefore,  Spark  Plasma 
Sintering  (SPS)  is  applied  to  a  mixture  of  NiTi-SE  and  Fe  powders  in  vacuum  at  various 
temperatures  for  a  period  of  5  minutes.  Since  University  of  Washington  did  not  have  an 
SPS  system,  when  this  project  was  started  in  2001,  the  SPS  located  at  Material 
Processing  Department,  Tohoku  University  was  accessed.  Later,  thanks  to  the  DURIP  AF 
grant  on  equipmental  proposal,  we  were  awarded  a  new  DURIP  grant  by  which  we 
purchased  a  new  SPS  machine  (Sumitomo  Coal  Mines,  Japan,  model  Dr.  Sinter  SPS- 
1020S).  Schematic  of  SPS  machine  is  shown  in  Figure  22.  The  interface  of  the  composite 
is  found  to  be  free  of  inter-metallic  compounds  which  are  known  to  the  make  the 
composite  brittle.  Profs.  Kang  and  Kawasaki  of  Tohoku  University,  Japan  helped  us 
process  the  FSMA  composites  used  in  this  early  part  of  the  work.  Though  the  current 
model  will  also  hold  true  for  laminate  composites  as  a  special  case,  particulate  FSMA 
composite  is  considered  due  to  the  superior  quality  of  composite  obtained  through  SPS 
processing.  Particulate  NiTi-Fe  composites  were  processed  by  SPS  at  600°C  in  vacuum 
for  5  minutes  and  50MPa  with  30,  50  and  70%  Fe  by  weight.  The  magnetization  (M)  - 
magnetic  field  (H)  curves  of  the  composite  specimens  are  then  measured  using  vibrating 
specimen  magnetometer  (VSM),  located  at  Tohoku  University,  are  shown  in  Figure  20. 
The  saturation  magnetization  (Ms)  taken  from  Fig.  23  are  compared  with  the  model 
results,  Fig.  18,  see  its  details  in  Chapter  3. 
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Figure  22.  Schematic  of  Spark  Plasma  Sintering  (SPS)  machine  used  to  process 
FSMA  particulate  composites. 


Figure  23.  M(emu/g )  -  H(Tesla)  curves  of  Fe-TiNi  composites  processed  by  SPS 

and  measured  by  VSM. 


As  far  as  the  magnetic  behavior  of  the  Fe-TiNi  composite  is  concerned.  Figure  23 
demonstrates  the  good  magnetic  behavior  that  is  predicted  by  the  new  model. 

We  tried  several  SPS  conditions  and  heat  treatments,  to  identify  the  best  SPS 
condition  (The  SPS  conditions  used  are,  600  C,  700  C,  and  800  C  all  under  50  MPa  for  5 
min,  in  vacuum). 
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Figure  24  (a)  and  (b)  show  the  DSC  results  of  Fe  (33.3  %wt  average  dia=100mm)  - 
TiNi  (average  dia=212pm)  with  these  three  different  SPS  temperatures  without  and  with 
aging  heat  treatment  of  320  C  for  30min. 


Aging  with  320C,  30  min. 


(b)  320  C,  30  min  heat  treated 

Figure  24.  Comparison  51  at%NiTi  (212  microns)  +  33.33%wtFe  (100  nm)  Composite 
processed  using  SPS  at  600, 700  and  800oC  for  5  min,  50  MPa  in  vacuum. 
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Figure  24  indicates  that  the  best  SPS  temperature  is  700C,  50  MPa,  5  min  in  vacuum 
conditions,  and  the  aging  increased  the  transformation  temperatures  shown  by  peaks. 
XRD  analysis  of  the  as-SPS  processed  with  700  and  900  C  were  performed  and  the 
results  are  shown  in  Figure  25  (a)  and  (b),  respectively. 


10  20  30  40  50  60  70  80  90 

20  (Cu  K  a)  /  degree 


(a)  51at%NiTi  (212pm)+33.33%wtFe  (74  pm)  composite,  SPS  at  700oC,  5  min,  50  MPa 


10  20  30  40  50  60  70  80  90 

29  (Cu  K  a)  /  degree 


(b)  51at%NiTi  (212  pm)+33.33%wtFe  (74  pm)  composite,  SPS  at  900oC,  5  min,  50  MPa 
Figure  25.  XRD  of  the  composite  NiTi  obtained  from  different  SPS  processes. 
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Figure  25  reveals  that  use  of  higher  temperature  for  SPS  resulted  in  excess  amount  of 
C  which  presumably  came  from  the  graphite  mold.  Therefore,  use  of  lower  temperature  is 
desired.  However,  lower  temperature  SPS  process  resulted  in  higher  porosity  in  the 
samples.  The  effects  of  SPS  temperature  on  the  density  of  the  samples  are  summarized  in 
Table  2.  From  the  results  of  Table  2,  DSC  data  of  Figure  24  and  XRD  data  of  Figure  25, 
use  of  700  °C  is  considered  to  be  best. 


Sample  ID 

SPS  . .  ..°C,  5min,  50  MPa 

Density 

measured* 

Theoretical 

density 

Porosity 

(%) 

600°C,51at%NiTi  (212  pm)+33.33%wt  Fe  (100  nm) 

6.66 

6.88 

3.2 

700°C,  51at%NiTi  (212  pm)+33.33%wt  Fe  (100  nm) 

6.65 

6.88 

3.2 

800°C,  51at%NiTi  (212  gm)+33.33%wt  Fe  (100  nm) 

6.65 

6.88 

3.2 

700°C,  51at%NiTi  (212  pm)+33.33%wt  Fe  (74  pm) 

6.56 

6.88 

4.6 

700°C,  51at%NiTi  (MRi)+33.33%wt  Fe  (100  nm) 

6.32 

6.88 

8.1 

900°C,  pure  51at%NiTi  (212  pm)  sample 

6.26 

6.40 

2.1 

600°C,  Dr.  Kang’s  sample 

6.07 

6.88 

11.8 

Table  2.  Density0  of  samples 

We  also  attempted  nano-sized  TiNi  powders  which  were  supplied  from  Materials 
Research  International  (MRi)  and  Argonide  Inc.  Figure  26  shows  the  DSC  results  of 
100%  TiNi  sample  processed  by  SPS  and  using  nano-TiNi  powders  (Ti  -51  at  %Ni  )  at 
850  °C,  5  min,  50  MPa,  followed  by  aging  at  320  °C  for  several  different  hours.  It  is 
obvious  from  Fig.  26  that  all  DSC  charts  did  not  exhibit  any  peaks,  indicating  no  phase 
transformation  during  heating  and  cooling  curves. 


*  Using  Archimedes  Principle 
pNiTi— >6.4  gm/cm3;  pFe— >7.86  gm/cm3 
0  in  gm/cm3 
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Figure  26.  DSC  chart  of  100%  TiNi  sample,  aging  at  320  °C  for  several  different  hours 


We  also  Processed  30wt%Fe  -  nano  TiNi  (10  nm  dia)  composite  by  SPS  at  800  °C, 
for  5  min,  under  50  MPa,  in  vacuum,  which  was  subjected  to  XRD  analysis.  The  result  of 
XRD  analysis  of  30wt%Fe  -  TiNi  composites  are  shown  in  Figure  27  as  open  symbol 
where  the  XRD  result  of  nano-TiNi  powders  (Ti  -51  at  %  Ni)  are  also  shown  as  a 
reference.  Figure  27  indicates  the  formation  of  excess  amount  of  TiNb  which  may  have 
suppressed  the  formation  of  TiNi  phase,  i.e.  B2  and  B19’.  Due  to  highly  reactive  nature 
of  the  nano-sized  TiNi  powders,  the  as-SPS  processed  Fe-TiNi  (MRi)  composites  did  not 
give  rise  to  any  SMA-SE  properties. 


Figure  27.  Comparison  of  XRD  results  of  TiNi  nano  powder  and  TiNi/Fe  composite 


Therefore,  we  focuse  only  on  use  of  TiNi  powders  whose  average  diameter  is  212  pm, 
supplied  by  Sumitomo  Metals,  Japan.  In  the  following,  all  TiNi  designation  referred  to 
this  larger  sized  TiNi  powder  (212  pm)  as  a  starting  material. 

The  cross  section  of  as-SPS  processed  30  %wt  Fe-TiNi  composite  is  taken  by  SEM 
with  EDAX  element  analyzer  at  UW-MSE  laboratory,  and  the  results  are  shown  in  Fig.28 
where  (a) ,  (b)  and  (c)  denote  the  cross  section  with  three  different  color  coded  elements, 
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Figure  28.  (b)  Larger  magnification  view  of  the  cross  section  of  Fe  (30  %wt)  -  TiNi  Composite 


200  jum 


Figure  28.  (a)  Cross  section  of  Fe  (30  %wt)  -  TiNi  Composite  where  particulate  phase  is  TiNi 

and  matrix  phase  is  Fe 
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Figure  28.  (c)  Line  profile  of  Fe  and  TiNi  phase  between  two  adjacent  TiNi  particulates. 


The  Fe  (30  %wt)  -  TiNi  (ave.  dia  =  212pm)  composites  that  are  processed  by  SPS 
with  700  °C,  5  min.  50  MPa  in  vacuum  are  subjected  to  compression  testing  at  different 
testing  temperatures.  The  results  of  compressive  stress-strain  (SS)  curves  are  shown  in 
Figure  29  where  (a),  (b),  (c)  correspond  to  the  SS  curves  tested  41°C,  -4°C,  -20°C, 
respectively.  The  SS  curves  tested  at  higher  temperature  show  superelastic  (SE)  behavior, 
while  the  SS  curves  at  lower  temperature  exhibit  shape  memory  effect  (SME)  behavior. 
This  is  due  to  the  fact  that  Af=  40  °C,  As=  -10  °C,  Ms=  -30  °C,  Mf=  -50  °C,  as  shown  in 
the  DSC  chart,  see  Figure  29  (d). 
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Chapter  5.  Processing  of  FSMA  Laminated  Composites  by  Plasma 
Etching/Cladding  Method 


Once  the  optimum  microstructure  of  a  FSMA  composite  is  identified  by  the  above 
modeling,  the  next  key  step  in  the  current  AFOSR  research  is  processing  of  such  a  FSMA 
composite.  During  the  past  several  years,  we  attempted  the  processing  of  FSMA 
laminated  composites  by  plasma-etching  and  cladding  method  where  three  kinds  of 
FSMA  composite  material  systems  are  focused  on: 

(1)  TiNi-FeCoV  with  A1  as  adhesive  layer 

(2)  TiNi-FeCoV  with  Ag  as  adhesive  layer 

(3)  CuAlMn-Fe  with  Mo  as  adhesive  layer 


Step  1 
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Figure  30.  Schematic  of  plasma-etching/cladding  process  (Taya,  2005). 


The  processing  of  a  FSMA  laminated  composite  by  the  plasma-etching/cladding  in 
vacuum  was  performed  by  the  facility  at  Toyo  Kohan  Co.  Ltd.,  Japan  as  this  company 
has  such  equipment.  The  concept  of  the  plasma-etching/cladding  is  shown  in  Figure  30. 
The  plasma-etching/cladding  has  been  successfully  applied  to  the  bonding  of  two  thin 
metal  plates  for  use  in  electronic  packaging  applications  (Taya,  2005).  However,  the 
order  of  the  thicknesses  of  ferromagnetic  plate  (FeCoV  in  the  first  attempt)  and  SMA-SE 
plate  (TiNi)  is  0.2-1  mm  range.  Therefore,  we  have  to  use  an  adhesive  thin  metal  foil  to 
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bond  these  metal  plates  after  heating  to  500  °C  range.  The  cross  section  of  the  above 
composite  plate  is  shown  in  Fig.  31  where  thin  aluminum  layer  was  used  as  adhesive 
material. 


Figure  31.  Cross  section  view  of'FSMA  composite  made  of  NiTi  (SME-SE)  and 
FeCoV (ferromagnetic  material)  with  aluminum  adhesive  layer,  processed  by  plasma-etching/cold 

rolling  process. 


It  is  noted  in  the  design  of  FSMA  composite  of  Fig.  31  that  use  of  FeCoV  is 
advantageous  as  a  soft  magnetic  material  because  its  saturation  magnetization  (Ms)  is  as 
high  as  2.2  T.  However,  to  make  FeCoV  converted  to  such  a  high  performance  soft 
ferromagnetic  grade,  we  have  to  apply  high  temperature  heat  treatment  in  H2 
environment  while  intermediate  heat  treatment  is  required  of  TiNi  to  convert  it  to  SE 
grade  in  air.  Thus,  we  have  to  make  a  good  compromise  in  selecting  a  combined  heat 
treatment  to  be  applied  to  FeCoV/TiNi  composite. 

In  order  to  improve  the  high  temperature  stability  of  TiNi/FeCoV  composite,  we 
applied  Ag  thin  adhesive  layer  instead  of  A1  for  bonding  of  TiNi  and  FeCoV  plates  as  Ag 
can  sustain  higher  temperature  than  A1  during  the  past  heat  treatments.  Fig.  32  shows  as- 
bonded  NiTi/Ag/FeCoV  composite  where  (a),  (b)  denote  the  overview,  cross  section  of 
the  composite  plate,  respectively  and  (c)  exhibits  the  EDX  analysis  of  NiTi/FeCoV  sheet 
composite.  DSC  measurements  and  B-H  curve  of  this  composite  after  it  was  heat-treated 
(800°C  30min  in  vacuum,  water  quenched)  are  shown  in  Fig.  33  (a)  and  (b).  Fig.  33 
indicates  that  the  NiTi/FeCoV  FSMA  composite  processed  by  the  plasma 
etching/diffusion  bonding  with  Ag  layer  exhibits  phase  transformation  under  temperature 
change  and  reasonably  good  magnetization  (M)  vs.  magnetic  field  (H)  relation. 
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There  are  no  inclusions,  compounds  at  interface. 


Figure  32.  TiNi/FeCoV  FSMA  composite  plate  where  adhesive  thin  Ag  layer  was  used  for  firm  bonding. 
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NiTi/FeCo_Nol_WQ  NiTi/FeCo_Nol_WQ 


Ms=-25°C,  Mf=-40°C,  As=-14°C,  Af=0°C 

(a)  (b) 

Figure  33.  (a)  DSC  data,  (b)  M-H  curve  of  TiNi/FeCoV  composite. 

Finally  we  attempt  to  process  CuAlMn-Fe  laminated  composites  by  plasma-etching 
cladding  method  where  Mo  was  used  as  adhesive  layer.  Figure  34  (a),  (b)  and  (c)  show 
as-bounded  composite,  the  cross  section  of  CuAlMn(bottom)-Mo(middle)-Fe(top),  at  low 
and  higher  magnifications,  respectively.  This  figure  indicates  strong  bounding  between 
them. 


Figure  34.  (a)  Appearance  of  CuAlMn/Fe  as-cladded  and  heat-treated 
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(b)  (c) 

Figure  34.  Cross  section  of  CuAlMn/Mo/Fe  composite  (b)  entire  thickness 
(c)  magnified  view  of  the  interfaced  area 
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Appendix  A:  Relation  Between  Bending  Moment  and  Curvature  for  Laminated 
Composites 

The  relation  between  the  normalized  bending  moment  and  the  normalized  curvature  of  the 
FMSA  composite  plate  is  classified  into  the  following  eight  patterns  as  shown  Figure  A.l. 


(a)  Case  1  (b)  Case  2  (c)  Case  3 

Figure  A.  1  Relation  between  Normalized  Bending  Moment  and  Normalized  curvature. 

Case  1  is  constructed  with  only  Pattern  1  (Figure  A.  1(a)) 

Case  2  is  constructed  with  pattern  1  and  2  for  the  loading,  and  pattern  1,  4,  5,  and  6  for  the 
unloading.  (Figure  A.  1(b)) 

Case  3  is  constructed  with  pattern  1,2  and  3  for  the  loading,  and  pattern  1,  4,  7,  and  8  for  the 
unloading.  (Figure  A.  1(c)) 

Equations  of  each  pattern  are  shown  as  follows. 

Pattern  1  (Case  1,2,3) 


where,  £  is  the  distance  of  the  neutral  axis. 
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Pattern  2  (Case  2, 3) 
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where,  £  is  the  distance  of  the  neutral  axis,  and  y,  is  the  position  for  a  =  a0 
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Pattern  3  (Case  3) 
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Pattern  4  (Case  2, 3) 


where,  £,  is  the  distance  of  the  neutral  axis,  and  y2  is  the  position  for  a  =  c r0  • 
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Pattern  5  (Case  2) 


where,  is  the  distance  of  the  neutral  axis. 


h  A, 


where,  is  the  distance  of  the  neutral  axis. 
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Pattern  7  (Case  3) 


where,  is  the  distance  of  the  neutral  axis. 
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Pattern  8  (Case  3) 
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where,  is  the  distance  of  the  neutral  axis. 
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USEFUL  RANGE 

The  useful  range  of  the  curvature  of  each  pattern  is  shown  as  follows. 
Case  1 

Pattern  1  (Loading  and  Unloading) 


Case  2 

Pattern  1  (Loading) 


Pattern  2  (Loading) 


Pattern  1  (Unloading) 
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Pattern  4  (Unloading) 


Pattern  5  (Unloading) 


Case  3 

Pattern  1  (Loading) 


Pattern  2  (Loading) 
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